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The renormalization algorithm based on regularization methods with two regula- 
tors is analyzed by means of explicit computations. We show in particular that 
regularization by higher covariant derivative terms can be complemented with di- 
mensional regularization to obtain a consistent renormalized 4-dimensional Yang- 
Mills theory at the one-loop level. This shows that hybrid regularization methods 
can be applied not only to finite theories, like e.g. Chern-Simons, but also to di- 
vergent theories. 



The renormalizability of a non-abelian gauge theory is not tied to the existence of a 
regularization method preserving its BRS symmetry, but to the more fundamental 
fact of the BRS symmetry not being anomalous [1]. Yet, having gauge invariant 
regularizations at our disposal is very useful for practical purposes. Sadly enough 
there are not many BRS invariant regularization methods available, dimensional 
regularization [2] and the method of higher covariant derivatives [3,4] being the 
two most prominent ones. 

Dimensional regularization is the most popular choice to regularize vector-like 
theories since the algebraic structure of these is not altered by a change in the 
dimension of spacetime. Ward identities hold almost automatically in this case. 
Unfortunately, when chiral objects like 75 matrices or Levi-Civita pseudotensors 
occur, an algebraically consistent dimensional regularization [5] is rather awkward, 
should it exist at all [6,7]. The reason for this is that the properties of chiral 
objects depend on the dimension of spacetime, which clashes somehow with the 
ideas behind dimensional regularization. 

The method of higher covariant derivatives [3,4] is a very interesting attempt 
to regularize non-abelian gauge theories while keeping the spacetime dimension at 
its physical value. Think for example of Yang-Mills theory. By adding suitable 
gauge invariant terms to the classical action -see eq. (4) below- one transforms 
the theory into a superrenormalizable gauge invariant theory, henceforth called 
higher covariant derivative theory. In this way partial regularization of Yang-Mills 
theory is achieved in a manifestly gauge invariant fashion. To regularize the few 
diagrams in the higher covariant derivative theory that are primitively divergent 
and at the same time preserve gauge invariance, one must supplement an additional 
gauge invariant regulator. Since these diagrams happen to be one-loop, finding a 
gauge invariant regulator that renders them finite should in principle be easier than 
addressing the regularization of the original Yang-Mills theory. The proposal in ref. 
[4] for this additional regulator is a certain explicitly gauge invariant Pauli-Villars 
regularization whose latest version is explained in ref. [8]. From now on, we will 
call higher covariant derivative Pauli-Villars regularization to this combination of 



higher covariant derivatives and Pauli-Villars regulators [8]. 

It has been recently shown [9] that the latter combination leads to a Yang- 
Mills theory with the wrong beta function. The reason being that there are non- 
local contributions to the vacuum polarization tensor that modify the non-local 
gauge-invariant part of the effective action. This no-go result renders higher co- 
variant derivative Pauli-Villars regularization rather useless as a general regular- 
ization method. Nonetheless, since the anomalous contributions come from the 
Pauli-Villars determinants themselves [9], it might be possible that some other 
regularization of the higher covariant derivative theory led to a sensible Yang-Mills 
theory. 

The purpose of this letter is to show that if we furnish the higher covari- 
ant derivative theory with dimensional regularization, one obtains a regularization 
method leading to a sensible 4-dimensional Yang-Mills theory at the one-loop level. 
We will call this regularization method higher covariant derivative dimensional reg- 
ularization and see that it is free of the problems that occur for higher covariant 
derivative Pauli-Villars regularization. 

In this state of affairs, one may rightfully ask why not using dimensional reg- 
ularization from the very beginning. In so doing we would save quite a lot of 
cumbersome algebra coming from the higher covariant derivative terms. In what 
follows we provide reasons to motivate the regularization method under consider- 
ation. Let us first consider perturbative three-dimensional Chern-Simons theory. 
The simplest higher covariant derivative term that one can add to the Chern- 
Simons action is F 2 /A, where A is a mass. This results into a very interesting 
theory on its own, namely topologically massive Yang-Mills theory, that when 
dimensionally regularized allows explicit two-loop computations. These computa- 
tions have been performed in ref. [7] and have shown the vanishing of the two-loop 
correction to the bare Chern-Simons parameter. They have also proved that topo- 
logically massive Yang-Mills theory is finite to all orders in perturbation theory, 
even though by power counting the theory is only superrenormalizable. With 



en eye on four- dimensional gauge theories involving 75 and the Levi-Civita sym- 
bol e^ u P a j it is therefore important to see whether the regularization method in 
question may become a general regularization procedure or it only works in three 
dimensions. Another reason is that if one shows that higher covariant derivative 
dimensional regularization is consistent, one then has an argument to not relin- 
quish the possibility of a complete 4-dimensional regularization based on higher 
covariant derivatives. Last, but not least, is that by carrying on with this reg- 
ularization method we shall learn how to renormalize theories whose divergences 
are parametrized by two regulators. This problem has not been tackled yet in a 
satisfactory way in the existing literature [8,10] even at one loop. 

Let us start by formulating higher covariant derivative dimensional regulariza- 
tion for 4-dimensional Yang-Mills theory. The classical Yang-Mills action for gauge 
group SU(N) in a covariant gauge dA a = b a is given by 

S = Sym + Sqf , (1) 



where 



Sys! = 1 jdScF^F^ 
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The action S is invariant under BRS transformations 



sA a = (D^c) a sb a = sc a = b a 



sc 



cabcJ>„c 



-gf aoc c°c c sJ« = sH a = 0, 



(3) 
the BRS operator s being nilpotent, s = 0. We have introduced the external 
fields J" and H a to keep track of the renormalization properties of the composite 
operators sA^ and sc a , respectively. To achieve regularization of the ill-defined 
Feynman diagrams corresponding to the action S in eq. (1), one proceeds in two 



steps [11]. First, one introduces higher covariant derivartive terms and replaces 
the classical gauge-fixed action S with the action 

S\ = SVm + S 2 + S gf , (4) 

where SVm is as i n ec i- (2) and S 2 and Sq F are given by 



S 2 = ^Jd*x(D 2 F fI1/ r(D 2 F^) a 
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Here a is the gauge fixing parameter and /(<9 2 /A 2 ) is given in momentum space 

by 



f( p -) = i + ^ 

J \A 2 ) A 4 ' 



so as to ensure locality and render all a-dependent contributions finite by power 
counting [9]. In our notation, a = corresponds to the Landau gauge and a — 1 
to the Feynman gauge. The action Sa in eq. (4) is invariant under the BRS 
transformations in eq. (3) and gives rise to a superrenormalizable theory that we 
shall call higher covariant derivative theory. The free gauge propagator for this 
theory has the form 



Gt(p) = 5 ab A ' 



•IWV" pV + A 4) 



(2 ^ A 4 



The ultraviolet behaviour of G ab {ji) is not strong enough as to make finite by 
power counting all 1PI Feynman diagrams. However, it is very easy to see that 
only one-loop 1PI diagrams contributing to the two, three and four-point Green 
functions of the gauge field are superficially divergent. In other words, the higher 
covariant derivative theory defined by S^ is superrenormalizable. To regularize 
it, we use dimensionally regularization and thus achieve full regularization of the 
original Yang-Mills theory while preserving BRS invariance. 



Having set a manifestly BRS invariant regularization method, however impor- 
tant this might be, is not the end of the story. Ultraviolet divergences have to be 
subtracted if a renormalized theory is to be obtained. In the case of gauge the- 
ories this means that given a regularization method, an explicitly BRS invariant 
subtraction procedure must be devised. Recall in this regard that dimensional regu- 
larization and ordinary abelian Pauli-Villars regularization lead to well established 
Bogoliubov-type subtraction algorithms [5,12]. To the best of our knowledge, no 
subtraction procedure has been shown to work for any higher covariant derivative 
regularization of non-abelian gauge theories. It would appear that setting such a 
procedure at one loop is a very simple task. After all, if a regularization method 
is BRS invariant, removing from the one-loop Green functions those contributions 
which are stricitly UV divergent yields finite renormalized Green functions that 
verify the BRS identities at one loop. This being true, one must not forget that is 
not enough. General renormalization theorems in Lorentz covariant perturbation 
theory [13] imply that renormalized 1PI Green functions at one loop are unique 
modulo a polynomial in their momenta. Hence, if a regularization process modifies 
the non-local structure of the theory in the large cut-off limit, serious dificulties will 
arise. This is precisely what happens with higher covariant derivative Pauli-Villars 
regularization. In this case, renormalization at one loop can be accomplished in a 
BRS invariant manner and yet non-local gauge invariant radiative corrections get 
modified so as to yield an unphysical beta function [9]. Let us see that everything 
goes smoothly when higher covariant derivative dimensional regularization is used. 

Let $ collectively denote the fields AjJ, c a , c a , b a , J® and H a , and let e 
denote the dimensional regulator defined by D = 4 + 2e , D being the spacetime 
dimension. The regularized one- loop effective action Ti(Q,e, A) obtained with 
the help of higher covariant derivative dimensional regularization depends on two 
regulators, e and A . We want to fully characterize its UV divergent behaviour. 
To do this, we first report the one-loop values for the three 1PI functions that 
will allow us to do so. The vacuum polarization tensor H~,(p, e, A), the ghost 
self-energy Q ab (p, e, A) and ghost vertex Vfi (k,p, e, A) as computed with higher 



covariant derivatives and dimensional regularization are given by 
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Here "vt" stands for contributions that vanish when one takes the sequential limit 
e ^ 0, A ^ oo, k is the dimensional regularization mass scale and p^ and k^ 
are external momenta. The symbols co(a) and ooq{q) denote constants which do 
not depend on k nor A. The function V^ n (k,p, a) is finite, has mass dimension 1 
and does not depend on any of the dimensionful parameters k and A. A number 
of comments regarding these expressions are in order. 

Comment 1. To obtain eqs. (7)-(9), we have first performed a Laurent expansion 
around e = 0, then dropped the terms that vanish as e — ► and finally discarded 
the terms that vanish as A — »■ oo. To do all this we have used the techniques devel- 
oped in ref. [7]. Notice that this way to proceed implicitly assumes that to obtain 
renormalized 1PI Green functions the regulators e and A should be removed in 
this order: first e —>■ and then A — ► oo. This prescription for removing the reg- 
ulators is rather natural for the regularization method we are considering and the 
only well-defined one for Chern-Simons theory [7,14]. Proceeding the other way 
around, i.e. taking first A — »■ oo and then e — > 0, is rather messy and would even- 
tually require renormalizing a theory at fixed but arbitrary e, something which is 
not clear how to achieve. 

Comment 2. The vacuum polarization tensor in eq. (7) has two types of UV 
singularities, a simple pole at e = and a logarithmic divergence in A . The pole 



comes from the fact that the higher covariant derivative theory has a divergent 
propagator at one loop. The singularity in A stems from the divergent structure 
of Yang-Mills theory. Note in contrast that the ghost self-energy and the ghost 
vertex do not show a singular behaviour at e = 0. This is due to the the fact that 
for finite A these two Green functions are convergent by power counting at D = 4. 
As regards other Green functions, it is not difficult to see using general results for 
one-loop dimensionally regularized integrals [15] and the techniques in ref. [7] that 
(i) only the three and four-point Green functions for the gauge field will become 

singular at e = 0, the singularity being a single pole, and 
(ii) that any other superficially divergent one-loop 1PI Green function will de- 
velop UV divergent singularities in In A. 

Comment 3. Using the identity 

In (p 2 /A 2 ) = In (p 2 /M 2 ) + In (M 2 /A 2 ) , (10) 

where M is a mass scale (later on to become the renormalization mass scale), we 
see that UV divergent contributions in the Green functions (7)- (9) are polynomial 
in the external momenta. Note that the non-local contributions explicitly displayed 
in eqs. (7)-(9) are the same as those obtained with only dimensional regularization. 
Furthermore, the techniques in ref. [7] ensure that also the non-local part in 
V^ n (k,p, a) in eq. (9) is the same as for dimensional regularization, a fact that 
will be used later. 

We now turn to the computation of the UV divergent contribution to the 
one-loop 1PI functional Ti(§,e, A). Recalling that the regularization method is 
explicitly BRS invariant and following standard techniques [16], we obtain that the 
one-loop effective action ri($, e, A) satisfies the BRS equation 

A A ,Dr 1 ($,e,A)=0, (11) 

where Aa d is the linearized Slavnov- Taylor operator for the D-dimensional coun- 



terpart of the action S\ in eq. (4). The operator Aa,d can be written as 

SSP 5 



fsD 6S 2 






where S§ is the D-dimensional version of the higher covariant derivative action 
in eq. (5) and 



A D = fd D x 



SS D 5 5S D 5 5S D S 5S D 5 5S D 



is the Slavnov- Taylor operator for the gauge-fixed Yang-Mills action in D dimen- 
sions S . In the limit e — » 0, A — ► oo we may express Ti(<&,e,A) as 

ri($, e, A) = rf v ($, e, A) + rf n ($) + "vt" , 

where again "vt" denotes contributions that vanish as e — > 0, A — ► oo . All 
finite contributions in this limit are collected in rf n ($), whereas r^ lv ($,e, A) 
only contains contributions which are divergent. The functional I^ 1 v ($,e, A) is a 
four-dimensional integrated local functional of the fields and its derivatives with 
mass dimension 4 and ghost number 0. Its coefficients are of the form aie -1 + 
ai In A + c^e" 1 In A, with ai, 02 and 03 constants. Using the methods in ref. [7], 
it can actually be shown that 03 = 0. 

Let us look at eq. (11). Since any one-loop 1PI diagram with at least one 
J® insertion is finite by power counting at D — 4, J-dependent contributions to 
Ti($, e, A) are finite as e — »■ 0. Recalling that all superficially divegent 1PI Green 
functions diverge as In A for A — *■ 00 and noting that the action S2 has an overall 
factor 1/A 4 , we have 

I lD 5SP <$ri 

hm hm lax -r-. — - = 0. 






Taking then the limit e -> 0, A ->■ 00 in eq. (11), it follows that rf v ($,e,A) 



satisfies 
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Ar*n$,e,A)=0 



where A is the Slavnov- Taylor operator for the gauge-fixed Yang-Mills action S 
in 4 dimensions. The most general solution to this equation has long been known 
[16]; it has the form 



-idiv 



rn$,e,A) = ciS Y M + AX, 



(12) 



where the functional X is given by 



X = / <fx 



c 2 {j;-d^)A a ^ + csH- 



a c a 



(13) 



and ci, C2 and C3 are coefficients. The coefficients c\, 02 and C3 can be obtained 
from the singular contributions to three suitable 1PI Green functions. Choosing for 
these the vacuum polarization tensor, the ghost self-energy and the ghost vertex 
in eqs. (7)-(9) and using the identity (10), we find that 
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Putting everything together, we have that r^ lv ( ( l ) , e, A) is given by eqs. (12)-(13) 
with coefficients as in eq. (14). 

Once we have characterized the UV divergent behaviour of ri(<&, e, A), we are 
ready to define a renormalized BRS invariant 1PI functional T r ^ n {^,g, M) that 
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describes Yang-Mills theory at the one-loop level. We define 

rr(^,g,M)= lim lim(l-T)ri($,e,A), (15) 

where 

TTi($,e, A) = (ci + /i) ^ M + A B 1 D , 

the functional X reads 

X D = jd D X [(C2 + / 2 ) (•/£ - ^C°) A a " + ( C3 + / 3 ) #°C ' 

and /i, /2 and /3 are arbitrary (e,A)-independent constants. Note that eqs. 
(12)-(14) guarantee that the limit in eq. (15) is finite. Taking into account our 
considerations above, it is very easy to realize that the equations 



A D TT l( *,e,A)=0 lim hm J dP x g ^_ T r l($ , e , A) = 



f 
hold, thus having 



lim lim A A ^rr($,^,M) = Arf n ($,^,M) 

A^oo e— >0 

and 

Ar\ en (<f>, g,M) =0. 

In other words, the renormalized one- loop functional defined in eq. (15) is BRS 
invariant. In fact, this is the most general 1PI functional that can be constructed 
from the regularized functional ri($,e, A) by performing local subtractions con- 
sistent with BRS invariance. This is due to the well-known fact [1,16] that the most 
general finite renormalization compatible with BRS invariance amounts to adding 
to the renormalized 1PI functional the quantity f\ S+A Jdzc [f2(J—dc) A+f^Hc\ , 
where j\, f'2 and fs are arbitrary finite constants. But this is already taken into 
account in the definition of the Bogoliubov-type operator T. 

11 



We now come to the computation of the one-loop beta function of the renor- 
malized theory. It is plain that the renomalized 1PI functional can not be obtained 
by multiplicative renormalizations of the coupling constant, the gauge parameter 
and the fields of the action S\ in eq. (4). The reason is that such a multiplicative 
renormalization will give rise to new divergences at e = since the higher covari- 
ant derivative term S2 will introduce divergent counterterms not associated with 
any primitively divergent diagram. Note that if we set a to zero and allow for a 
renormalization of A, or equivalently multiply S2 by a bare parameter Ao, our 
renormalized 1PI functional can be obtained by multiplycative renormalization of 
a regularized bare theory having S\ as classical action [10]. Strictly speaking, 
this renormalization procedure demands three renormalization conditions to de- 
fine Yang-Mills theory, for one has to define what is meant by renormalized A, or 
equivalently renormalized Ao . It is well-known [1] however that the construction of 
Yang-Mills theory along the lines of renormalized perturbation theory [13] requires 
only two renormalization conditions, as is the case with dimensional regularization. 
Hence, if we did renormalize A, we would have to understand Yang-Mills theory 
as the large-A limit of a renormalized theory having S\ as classical action. In this 
paper we are not concerned with this approach -it fails altogether for q^O [10]. 

Having no multiplicative connection between the renormalized theory and a 
regularized bare theory, standard text book techniques can not be used to com- 
pute beta functions and anomalous dimensions. Yet this should pose no problem, 
since the knowledge of the renormalized Green functions is all that is needed to 
fully characterize the theory. Let us recall in this regard that for the rigorously 
established BPHZ subtraction procedure there is no bare theory [1], the renormal- 
ized theory is obtained directly and everything, if liable of a perturbative calcu- 
lation, can be computed from renormalized quantities. Since the beta function 
and the anomalous dimensions are the coefficients in the renormalization group 
equations for the renormalized Green functions, they can be determined from the 
latter equations. In our case the renormalization group equation for a 1PI Green 
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function Gft(p e ,g,a, M) reads 
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G R (p e ,g,a,M) = 0. (16) 



The renormalized vacuum polarization tensor, ghost self-energy and ghost vertex 
in the minimal scheme fa = fa — fa — are given by 
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(17) 

Substituting these expressions in eq. (16) and expanding in powers of g the 
functions (3(g), 7$(g) and 5(g) , one obtains a linear system to be verified by the 
coefficients of such series. Solving this system leads to 
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These results are in complete agreement with the beta function and the anoma- 
lous dimensions computed within any consistent mass independent renormalization 
scheme known as yet, say the MS scheme for dimensional regularization. 

Let us make two final remarks. First, note that the non-local part of the 
renormalized Green functions in eqs. (17) are the same as for any renormalization 



• Note that the renormalized 1PI functional T\ cn (&,g, M) is defined by a Bogoliubov-type 
subtraction from a regularized 1PI functional compatible with a quantum action principle. 
This ensures that the renormalization group equation holds true for the renormalized Green 
functions. 
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scheme based on dimensional regularization. This, together with the fact that the 
renormalized 1PI functional can be constructed from the three Green functions in 
eqs. (17), implies that the renormalized effective action r^ en (<^>, c/, M) we have con- 
structed here and the one-loop dimensionally renrormalized 1PI functional differ by 
finite renormalizations of the fields and the parameteres in the theory. Our second 
remark concerns the existence to all orders in perturbation theory of a consistent 
BRS invariant renormalization algorithm based on higher covariant derivative di- 
mensional regularization. Although this problem lies well outside the scope of this 
paper, let us say that the subtraction procedure should take care of the singular- 
ities at e = first and then remove the divergent large-A behaviour. The BRS 
invariance of the algorithm should come from BRS invariance of the regularization 
method and the fact the the purely divergent UV contributions are BRS invariant 
on their own. Although we do not have a rigorous proof of these statements, we 
have partially checked them for two-loop diagrams with overlapping divergences. 
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